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ABSTRACT 

We have extended the theoretical considerations of Scudder and Olbert 

(1979) (hereafter paper I) to show from the microscopic characteristics of 

the Coulomb cross section that there are three natural subpopulations for 

plasma electrons: the subthermals with local kinetic energy E < kT^; the 

transthermals with kT _ < E < 7 kT and the' extrathermals E > 7 kT . We 

C 0 ' ■ c 

present experimental support from three experimental groups on three 
different spacecraft in the interplanetary medium over a radial range for 
the five interrelations projected in paper I between solar wind electron 
properties and changes in the interplanetary medium: 1) subthermals 

respond primarily to local changes (compressions and rarefactions) in 
stream dynamics; 2) the extrathermal fraction of the ambient electron 
density should be anti -correlated with the asymptotic bulk speed; 3) the 
extrathermal "temperature" should be anti-correlated with the local wind 
speed at 1 AU; 4) the heat flux carried by electrons should be 
anti-correlated with the local bulk speed; and 5) the extrathermal differ- 
ential "temperature" should be nearly independent of radius within 1 AU. 

From first principles and the spatial inhomogeneity of the plasma we 
show that the velocity dependence of Coulomb collisions in the solar wind 
plasma produces a bifurcation in the solar wind electron distribution 
function at a transition energy , E* . This energy is theoretically shown to 

scale with the local thermal (j' core) temperature as E*(r) (Z* 7 kT (r). 

c 

This scaling is observationally supported over the radial range from 0.45 
to 0,9 AU by Mariner 10 data and IMP data acquired at 1 AU. The extra- 
thermals, defined on the basis of Coulomb collisions, is synonymous with 



3 


the subpopulation previously labeled in the literature as the "halo” or 
"hot" component. If the transition energy should be required to equal the 
polarization potential energy, e'^(r) , for all radii beyond the observer, 

(as motivated by Mariner 10 data) the thermal electrons should obey a 
polytrope law with index Y = 7/6 . In this circumstance the polarization 
potential is equal to the specific enthalpy: e$(r) = [y/(y- 1)] kT^(r) . 

This relation probably does not obtain in the corona, inside r < 10-20 
(cf. Paper I). In the asymptotic spherically symmetric solar wind the 
inverse power law index for radial variation required for the thermal 
electrons should be a = 1/3 which is consistent with the recent in situ 
determinations between 0.45 and ^3 AU. We thus provide the first 
self-consistent argument for associating E*(r) with e$(r). This 
theoretical asymptotic (r -► “) thermal electron temperature variation is 
between the conduction dominated (2/7) and the inviscid solution (2/5). In 
the proximity of the ecliptic the ratio of the thermal electron Coulomb 
mean free path to scale length is shown to decrease with increasing radial 
distance. By contrast, over the solar poles the same asymptotic tempera- 
ture and density variations imply that the. Coulomb theTmal mean free path 
will increase with increasing radial distance. 

A straightforward extension of our model to slow time dependent 
situations shows that the fundamental scaling of E* with temperature and, 
by inference, the polytrope law for thermal electrons is theoretically 
maintained within compression and rarefaction regions. An illustrative 
calculation is given of the extrathermal density's response in stream 
compressions near I AU. These calculations show that the extrathermal 
electrons density can either increase or decrease depending on whether the 
core electrons compress more nearly as an isothermal gas or an adiabatic 
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one. On the basis of these calculations and a reexamination of published 
data we suggest that the extrathermal density is probably enhanced within 
stream compression regions. This implies (as observed) that (a) the core 
electrons should compress more nearly as an isothermal (y = 1) rather than 
as an adiabatic (y = 5/3) gas, and (b) a very large electron heat flux- 
should be observed leaving the compression region along the magnetic field. 

We thus conclude that the global and local Coulomb effects discussed in 
paper I and elaborated here are essential aspects of the solar wind plasma 
as it is observed. The experimental support which the predictions of the 
theory in paper I enjoy, strongly suggests that the further development of 
that theory is warranted with the ultimate objective of understanding the 
impact that global heat transport has on the physics of coronal expansions . 
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INTRODUCTION 

There is, even now, no self-consistent theoretical understanding of how 
the thermal and/or mechanical (wave) energy of the lower corona supports 
the spectrum of solar wind states observed in the interplanetary medium. 

The usual theoretical simplifications for wave propagation (WKB) and heat 
conduction (small ratio, K, of the mean free path' to the scale height known 
as the Knudsen number) cannot be defended at crucial phases of the 
acceleration of the incipient wind. 

To our knowledge there exists no fully self-consistent theory for heat 
transport in a Knudsen regime (K = 1) plasma such as is found immediately 
above coronal holes from which the (high speed?) solar wind is thought to 
emanate. Various ad hoc localized collisional models have been proposed 
(cf. Rawls et al. . 1975) to allow generalizations (for K “ 1) of the local 
Spitzer (1962)-Braginskii (1965) result Q = - kVT, often used in solar wind 
models, which only rigorously applies for K << 1 in nearly homogeneous' 
systems. We have recently suggested (Scudder and Olbert, 1979, Paper I) 
that the transport physics for solar wind electrons (who carry most of the 
heat flow) in the Knudsen regime should have a global character alien to 
these localized models; this new theory has suggested several inter- 
relationships between experimental electron parameters in an inhomogeneous, 
fully ionized Knudsen regime plasma that should be implicit in solar wind 
electron data. 

The considerable, additional effort required to make this theory 
mathematically self-consistent (to allow accurate -theoretical models of 
thermally driven coronal winds) is unwarranted unless these inter- 
relationships and other predictions of the theory are supported by the ever 
growing data base of solar wind electron measurements. This comparison 
represents an extensive test of the theory since these observations now 
span the radial domain of 0.3 ~ 5 AU and represent more than a decade of 
monitoring at the orbit of earth. If, however, the theory in its present 
form predicts more than it assumes, further development of this novel, 
global theory of Knudsen regime heat transport would appear fruitful. 

This paper develops the predictions of the theory and compares them 
with the relevant solar wind electron observables. It represents a review 
of the salient solar wind electron morphology as organized by the 
theoretical ideas of paper I. Our objective is to evaluate the merit of 
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that approach to the transport physics of plasma electrons in the Knudsen 
regime. Our conclusion is that the theory, although lacking rigor, 
predicts and organizes a large number of reported electron observations; by 
virtue of this extensive experimental support the theory warrants further 
development . 

In the first section of this paper we develop from first principles the 
idea of subthermal (E < -kT^) , thermal (E = kT^) , transthermal (1 < E/kT^ < 
7), and e'xtrathermal (E > 7 kT^ = E*) electrons as determined by the 
dynamics of Coulomb collisions; (k is Boltzmann's constant and is the 
temperature) . The preceding definitions are based on the assumption that 
the subthermal electrons are Maxwellian. Statistical arguments identify 
the subthermal electrons, E < kT , as the electron subpopulation that has a 
greater than three sigma chance of having a Coulomb 90 cumulative 
scattering while traversing a (rms) thermal speed particle's free path. 

This subthermal portion of f^ contains the dominant fraction of the so 
called "core” density. The transthermal subclass kT < E < 7 kT does not 
have a clear local or global character. This energy regime is statisti- 
cally more local at kT^ and most global toward' 7 kT^. We demonstrate that 
the extrathermal subclass of electrons are those locally observed electrons 
which as a class have a greater than 95 % chance of having had their last 
•Coulomb collision more than one (rms) thermal electron's free path away 
from the observer. With experimental data taken at a variety of radial 
distances within and at 1 AU, we show that the breakpoint, or "transition” 
energy, E», in the solar wind electron distribution function scales as 
E»(r) - 7 kT^(r). This scaling is that expected for the theoretical lower, 
bound of the Coulomb defined extrathermal population. We thus identify the 
theoretical Coulomb extrathermal population with the previously described 
"halo” or "hot" component f Feldman et al. (1975)), and illustrate those 
aspects of Coulomb collisions which allow this subpopulation to be 
discernible in the solar wind electron velocity distribution function. In 
section 2 we discuss the relation of this transition energy E* to the 
magnitude of the electrostatic polarization potential and develop the 
evidence for the polytrope behavior of thermal electrons with polytrope 
index y - 7/6, 

These theoretical considerations of the properties of Coulomb 
collisions in a plasma which allow the "halo-core" dichotomy permit more 
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precise statements of the expected correlations in steady state made in 
paper I which are developed in section 3. For the convenience of the 
reader, we restate them here: 

(1) In contrast with extrathermal electrons, subthermal solar wind 
electrons should behave most nearly as a "classical gas". We define a 
classical gas as an ordinary gas with finite mean free path (mfp) effects 
(e.g., heat conduction) in contradistinction to an ideal gas where mfp = 0 
(and conduction vanishes). A classical gas will compressi-vely heat and 
expansively cool to the extent that heat transport allows. 

(2) The solar wind extrathermal fraction of the electron density should 
be anti-correlated with the local bulk speed (if it is nearly the 
asymptotic wind speed). 

(3) The extrathermal electrons should generally form a subpopulation of 
the electron component of the solar wind plasma characterized by a 
"differential temperature", T , distinguishable from T . The former 
temperature is expected to be anti-correlated with the local bulk speed at 

1 AU. The ratio 7^^/T^ is determined by the ratio of the average 
collisional depths to the corona of each subpopulation. 

(4) The heat flux carried by electrons should be anti-correlated with 
bulk speed. 

(5) The differential temperature of extrathermals should be nearly 
independent of radius in the inner heliosphere (r < 1 AU) . Section 4 
contains the experimental support from new data (Mariner 10), and from 
those previously published (IMP, Helios) for the above stated relations. 

In the final section we discuss a simple extension of our steady state 
model to the case of slowly evolving stream patterns which is the general 
situation at stream-stream boundaries. We calculate there the response of 
the halo density within compression regions and illustrate its sensitivity 
(through the induced gradients of the core electron properties) to the 
local modifications of the polarization potential. 

general features of the observed and theoretical electron VELOCITY 

DISTRIBUTION FUNCTION 

The solar wind electron velocity distribution function is generally 
observed at all radial distances to be nearly Maxwellian in shape at. low 
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energies, ("the so-called "core”) with a smooth monotonic transition to a 
hotter, sparser, but nearly isotropic suprathermal population, (the 
so-called "halo”). To illustrate this point four representative exfimples 
of cuts through the velocity distribution at different radial locations are 
displayed in Figure 1. Velocity distribution functions with these general 
shapes were theoretically calculated in paper I. The generality of those 
shaped solutions is not obvious , given that radial profiles had to be 
assumed in the "bootstrap" calculations presented there. By very general 
arguments we shall show in this section that such an electron velocity 
distribution function is almost unavoidable for a plasma locally in the 
Knudsen regime that is part of an inhomogeneous plasma system such as is 
found within the heliospheric cavity. In so doing we also indicate 
additional theoretical predictions of the expected scaling of the 
transition energy between thermal and suprathermal subpopulations. We 
shall contrast the theory with experimental data and show that the expected 
scaling is indeed observed. 

We suggested in paper I that the extrathermal subpopulation of solar 
wind electrons were last members of collisional distributions much nearer 
to the corona than the corresponding place for the thermal electrons. 
Qualitatively,, the situation is analogous to the situation in radiative 
transfer, where the same observer at different energies (frequencies) is at 
different collisional (optical) depths to the "source" of the phase density 
(emission) . 

"Survival" Probabilities 

The Coulomb collisional momentum transfer scattering frequency may be 
written as 

V . ‘•e- ‘•pi 

fg, fp) 

where u = w/w„ (with w being the proper frame speed and w„ = 

f 3 kT /m" ) , and where the collision frequency, v , and integral (u) are 
as defined in paper I. The inverse cubic dependence of v on "test" 

.particle proper frame ^peed is the well known two-body Coulomb result. The 
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dimensionless integral , accounts for the thermal motions of the electron 
and ion "targets" with which a "test" electron interacts. The speed 
dependence of tSl (u) is of some interest in defining natural subpopulations 
of electrons with regards to Coulomb collisions. This functional depend- 
ence is illustrated in Figure 2 for Maxwellian electron and ion distri- 
bution functions of nearly equal temperatures (T /T = ra /m ’ 0.93) 

G P • 6 P 

for the range of "test" electron speeds much less than the ion "target’s" 
(rms) thermal speed, w._ , to much larger than the electron "target's" 

9 P 

(rms) thermal speed, w,j, 

There are two contributions to (u) : 1) due to the "test" electron 

scattering off of target ions and 2 ) due to "test" electrons scattering off 
of target electrons. The cross hatched portion of Figure 2 indicates the 
electron-electron scattering contribution to^(u). The remaining distance 
below the curve results from electron-ion scattering. From this graph it 
is clear at a given speed anywhere on this diagram that e-p scattering, 
which is an elastic process, is the dominant contribution to ^ (u). Almost 


all solar wind electrons are suprathermal insofar as ion most probable 
speeds are concerned (i.e., w/w^ p» 1 ), which implies that the 
electron-proton portion of^(u) is :£ 1. However, as the stippled ribbon 
in this figure indicates, typical solar wind electron observations straddle 


the similar thermal to suprathermal transition of electrons scattering off 
of background electrons. For electrons of speeds w » j^2w„ (u) 

asymptotically approaches 3/2. Physically this arises, from the fact that 
the momentum transfer is most efficient for particles of comparable speed 
and there are at these speeds very few additional "target" electrons 
available. 


Using the probability concepts of paper I we desire to calculate as a 
function of proper frame speed the probability, P 1 • 

sca u uor* 

p - 1 - p 

scatter survival 

that a "test" electron will be scattered while traversing a slab of 
thickness AZ equal to the local Coulomb free path, computed for the 

(rms) thermal speed of the electron, g* We thus obtain 

^survival = Wp' {-[J (u)/Jl (1 )]/u^) 

which is plotted in Figure 3 . The solid curve indicates the probability of 

survival (right axis), [of being scattered (left axis)], as a function of 
scaled proper frame speed u = w/w^^^. The cross-hatched region indicates 
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the typical regime of solar wind plasma electron measurements at 1 AU. 
"Test” electrons with proper frame speeds greater than 
w» = u* w„ u* = 2.16, Off), 

will transit unaffected by Coulomb collisions with a greater than 3^ 

probability. Nevertheless there is a finite (< 5 %) probability of local 
interaction with the ambient "targets". We have called this statistically 
defined subpopulation, the "extrathermal" electrons, since their 
collisional antecedents are beyond those of the thermal population. 


t|? ”t, 


„ < w < u* w^, „ are the 

o 1 y € 


The electrons of intervening speeds 
"transthermals" and have intermediate chances of scattering (survival) as 
may be determined from Figure 2. 


Electrons with speeds 




have a greater than 3a chance of 


being scattered while traversing one These particles have proper 

frame energies less than kT and we have called them the "subthermals" . 

We now address the question of how the preceding calculations allow for 
a bifurcation of the electron distribution function. Thermodynamic equi- 
librium is a situation where extrathermal electron departures and arrivals 
to the observer are in detailed balance; in the strict sense of the word 
thermodynamic equilibrium can obtain only in a homogeneous medium. In this 
circumstance the shape of the distribution function is time independent 
even though individual extrathermal electrons are constantly moving 
quasi-ballistically further than the thermal electron's mean free path, 
Xmfpj between collisions. In contrast with this "homogeneous" thermo- 
dynamic equilibrium, the solar wind processes take place in a radically 
inhomogeneous plasma. As shown in Figures (5,6,7) of paper I there are 
strong gradients in the density and the smooth forces implied in the solar 
wind expansion which affect the evolution of the distribution function even 
in portions of the solar cavity where gradients are weak! Most astro- 
physical plasmas are inhomogeneous ultimately as a result of gravitational 
fields; this fact plus the Coulomb "window" above u* permits more dense 
extrathermals in one proximity to "leak" through to the local medium and 
show up as an overpopulated extrathermal population relative to the local 
collisional Maxwellian which has (by inhomogeneity) a different density 
and/or temperature. 
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Because of the properties of Coulomb collisions in a plasma and the 
inhomogeneous nature of the solar wind plasma, we therefore expect a 
noticeable change in the shape of the electron distribution function in the 
vicinity of the proper frame kinetic "transition” or breakpoint energy 

®*Iheory 


E* 


Theory 


^ (u»)^ kT^ s 7 kT^. 


( 1 ) 


Because the preceding definition has been made on the basis of the 
statistical probability of Coulomb interactions we should expect a 
dispersion of the experimental determination of u* . 

A transition energy to the breakpoint energy may be 

determined directly from the data. It has been customary to fit the 
thermal electrons and suprathermal electrons with two separate (bi-) 
Maxwellian profiles . ' This is empirically motivated by the sharp transition 
of slopes of the Inf^ vs. E plot between tvra straight segments which 
approximate, respectively, the sub /trans thermals and extrathermal electrons 
in our resent terminology. 

We now offer quantitative support from in situ measurements of the 
transition energy of solar wind electrons. 


Mariner 10. CMVM '7^1 

We have determined the transition energy from the data acquired by the 
FLasma Science feperiment (PSE) on Mariner 10 between 0.45 and 0.9 AU. 

This breakpoint energy has been operationally determined as that proper 
frame energy where the core and halo Maxwellian model profiles intersect. 
The core and halo fits were made (cf., Ogilvie and Scudder, 1978) to data 
on either side of the break in the observed distribution function; as such 
E*jjata intersection- of the extrapolated fit profiles. Nevertheless, 

E*jjata quantifies a region of relatively strong curvature of Inf vs. E, 
between thermal and suprathermal electrons. 

In Figure 4 we illustrate, in histogram form, all of the hourly 
averaged scaled breakpoint energies E*j^(r)/kT^(r) where T^(r) is the core 
temperature at radial distance r obtained by the Mariner 10 PSE. The 
functional scaling by kT is that suggested above by the Coulomb theory of 
the transition energy. The binned quantity is dimensionless, allowing data 
with disparate temperatures at various radial positions to be analyzed 
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together. The binning has also been done for disjoint radial subintervals 
with results similar to those shown in Figure. 4. We thus suggest that the 
superposition of observations from different radial positions gives a 
statistical description unbiased by radial variations. In particular the 
width is what an observer should see with comparable statistics at a fixed 
radius. The principle features of this histogram are its modal and average 
values : 




= 7.37 
modal 




kT^(r) 


= 7.87, 
ave 


its narrowness and its extreme values 

p. p ^ c 

5.5 .< <. 10.5. 

kT^(r) 

From these values we infer a spread of scaled speeds 

1.91 < u» < 2.64 

about most probable and average values of 


MVM 


=2.21 and u» 


modal 


MVM 


= 2.29. 


average 


The above numerical relations were determined between 0.45 and 0.9 AU. The 
modal value has been indicated by the arrow in Figure 2. The entire range 
'^*MVM given by the vertical stippled zone. This experimentally 
determined, scaled break speed nearly equals the theoretical value (2.16) 
discussed above where there is a three sigma (95^) chance of transiting the 
(rms) thermal electron's free path at the observer unaffected by Coulomb 
pitch angle effects . It is interesting to note that the experimental range 
of u* for all Mariner 10 data indicated by the vertical stippled zone in 
Figure 3 is very nearly centered (logarithmically) within the theoretical 
2o and 4a criteria of definition of the extrathermal subpopulation 
(horizontal, diagonally shaded, zone in Figure 3)- We interpret this break 
in the slope of Inf vs. E in the vicinity of E « 7 kT as the lower limit 

e .Q 

of the extrathermal component defined above. 
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^*Data-^ 

In order to shbvr that other experimental data at another spatial 
location give similar results consistent with our modal and intrinsic 
extreme values of reproduce in Figure 5 the correlation data 

of En vs. T from Feldman et al. (1975). These authors describe (p. 4192) 
how Eg was operationally determined : "The boundary between halo and core 

components ... is that energy above which the count rates first exceed the 
count rate predicted by the best fit function f^ [fit to the lower 
energies] by more than 1.5 standard deviations." ..."...Eg is the break 
point energy (placed midway between that energy level determined above and 
the next [sic] lower) averaged over the four energy sweeps (covering a 90° 


sector) which bracket the heat flux direction Q ." "The error in deter- 

e 

mining Eg is about +1/4 the energy level spacing. This is variable, since 
the IMP levels are logarithmically spaced, but for E„ 60V, o„t, “ 5V." 

D i^D 

The best linear model of the form E„ = A + B kT for these data has 

D C 

been overlayed in Figure 5. This line has the equation 

E„(eV) = -17.5eV + 7.25. kT^(eV) 

From this analysis the optimal variation of E„ vs. T is 

D C 


+2.95 

B = dE„/d(kT„) = 7.25 

D O 

-1,75_ 

with an intercept value of 


A . - 17.5 lll’l (eV) 


The systematic differences between Eg and E* determined via the procedure 
quoted above are of three types: 

(i) the difference between the break definitions of 
Feldman ^ al.. (1975) near and the break point 
EgA* which is the average of Eg over all viewing 
angles , 

(ii) the quantization of the energy intervals between 
data points in the breakpoiint energy interval, 
and 

(iii) the IMP breakpoint energy was determined in the 
spacecraft frame rather than the proper frame. 
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The effects of (i) are of the order [+3.1 ± 193 eV (cf., Table I Feldman ^ 
al . (1975) with Eg < . The effects of (ii) are of the order of one 

half the average energy interval of quantization present in the graph of 
Figure 5 — this systematic effect is +10 eV. At 400 km/sec and 
j'60 eV effect (iii) is a 13/6 uncertainty. Thus the fit intercept values of 
A are consistent with zero, when systematic effects are considered. The 
systematic uncertainty of the intercept is 13 eV. Thus Eg* is consistent 
with being strictly proportional to T^ with no offset. 

The extreme values of the IMP slope determination are consistent with 
those determined from the Mariner 10 PSE data. This additional 
determination gives 


IMP 


2.19 


1 AU 


+0.41 

-0.27 


which is certainly consistent with the theoretical 3a value of 2.16 


+0.64 
- 0.26 ’ 


where the theoretical dispersion arises from the 2 and 4a definition of 
E*Theory (sQuation 1). The variance of the slope of the IMP data is 
consistent with the Mariner 10 widths being intrinsic . Both data sets over 
a range of radial locations are consistent with the predicted spatial 
invariance of which is predicated-on the physics of Coulomb 

collisions . 

Thus, the above results suggest that the linear relationship between 
E*^ , and T is a consequence of Coulomb collisional dynamics. By 
contrast Feldman ^ (1975)' p. 4184 suggest that this 'relationship (1) 

"may be no more than the expected kinematic relation between Eg and (see 
Appendix A)", (where these authors show that given two Maxwellian fits to 
the core and halo and if T^/T^ = const then Eg kT^) , or (2) may arise 
because e* estimated from the electron energy equation is dominated by the 
enthalpy term. 


2. THE TRANSITION ENERGY -GLOBAL AND LOCAL CONDITIONS RECONCILED 

Until the calculations of paper 1 based on Coulomb collisions it had 
been thought that only (collisionless) exospheric theory could predict the 
change in the spectral shape of the electron velocity distribution function 
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at the energy we have called the transition energy. The "natural" expla- 
nation of this separatrix relied on the different expected behavior of two 
classes of collisionless electron orbits: 1) those that are energetically 

bound in the electrostatic polarization potential (the "core") and 2) those 
that are free in this potential, (the "halo"), (cf. Jockers (1970, Schultz 
and Eviatar (1972), and Perkins (1973)). Because the collisionless 
(Vlasov) exospheric problem does not provide an unambiguous solution (for 
example anv function of the constants of the motion is a solution) , it must 
be pointed out that the distribution functions selected by these authors 
were strongly influenced by the available experimental distribution 
functions. Given the theoretical and empirical discussion of the previous 
section including Coulomb collisions and the successful calculation of 
steady state distribution functions in paper 1 that closely resemble real 
data, it is fruitful to reopen the discussion of a transition energy which 
depends only on local properties (from a collisional point of view) with 
the transition (or breakpoint) energy which is identified in collisionless 
exospheric theory with the polarization potential energy. This potential 
depends in steady state on a global distribution of electron pressure, 
density and other variables along the entire tube of force. We shall 
develop the necessary and sufficient conditions that these two interpre- 
tations are compatible in steady state. This matching condition requires 
in steady state that the thermal electrons obey a polytrope law with 
specified polytropic index. The polytrope "equation of state" offers the 
possibility of a "local” description for the thermal electrons. The 
predicted value of the polytropic index predicts more nearly isothermal 
than adiabatic behavior as appropriate for thermal electrons. The theory 
makes a definite prediction of the radial variation of thermal electrons 
when the density variation is known. The consequences of this matching 
condition are shown to be. remarkably accurate by direct comparision with 
relevant experimental data. 

From the discussion of section 1 it is unclear what theoretical 
relationship, if any, electrostatic polarization 

potential $. Nevertheless or its relatives such as Eg, are often 

taken to be the polarization potential energy barrier to infinity. It has 
been suggested^ for example, that Eg responds coherently with estimates of 
the polarization potential inferred from the electron energy equation when 
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the solar wind flow is nearly steady and spherically symmetric. .. [see 
Feldman et al. (1975)]. 

In steady state the potential energy barrier to infinity is the work 
required to move an electron along any given magnetic field line to 
infinity; 


e($(~) -$(r )) - X Z e E • dn = Z e E. . dS,, (2a) 
P(r^)- / - II 

where P(r^) is the observation point and ^ indicates the line integral 
along the appropriate £ line. The last two integrals of (2a) are equal as 
a result of the path independence. If we ignore the additional electric 
field in the proper frame due to the ohmic term, propagating waves, and 
deviations from charge neutrality in the solar wind beyond 1 AU we have 
from the generalized Ohm’s law [see for example Rossi and Olbert (1970)]: 

Z-eE^VPg/n^, (2b) 

where P is the electron partial pressure, n„ the electron number 

- 10 ® 

density, where Z_ is -1, and e 4.8 x 10 esu. Since the electron 
pressure P^ is dominated by the thermal partial pressure, we may write the 
potential as 


eO(ro) = kT^(r^) 6(r^), 

where 


(3a) 


S(r^) 





7(ln n^)'d£ 


(3b) 


where we have assigned the zero of potential at infinity. 

The theoretical, polarization potential *(i^q) scales directly with the 
local temperature provided 5(r^) is a weak function of the local 
temperature. If this situation obtains, then 


6(ro) 


d(e$) 

d(kT^) 


r 


0 
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should empirically be a constant with respect to variations at a fixed 
radial position. 

From experimental data there is strong support for this idea. By 
combining the information in 1 AU data given by Feldman ^ al. ,(1975) 
(Figures 2 and 18) and noting that in these author's notation is e$ in 
ours, we see that 

« . , , I = V ^ 

experimental | 

1AD 


which is apparently not a function of local temperature, T . We have used 
the fact as argued by those authors that (E,p, Eg) and (Eg, kT^) 
linearly correlated. It is by no means clear that equation (3-b) 
transparently implies that 6 ( 1 ^^) is nearly independent of Tq(i^q)? however, 
this appears to be empirically substantiated — at least on a statistical 
basis. It is therefore clear from microscopic Coulomb properties and the 
mathematical structure of the polarization potential that at 1 AU in steady 
state E* and e$ should be proportial because of their common parametric 
dependence on the local temperature, T^. 

Mariner 10 data indicate that on average (r)/kX^(r) is 

independent of radius between Mercury an d Ea rth. It is therefore possible . 
given the evidence at 1 AU, that the Mariner 10 data imply that 


e$(r) 

kT^(r) 


= 5(r) 


is independent of radius (so long as the conditions outside of observer are 
similar to those outside of the 1 AU observer.) 

Pursuing the consequences of this hypothesis in mathematical terms, we 
should have d 5 (r)/dr “ 0. This condition is only satisfied if the core 
obeys a polytrope relation for all radii beyond that of the observer; 

P/Pq = (n/n^)^; Y = 6/(6-1)5 T^,(r)/T^(r^) = (n/n^) > *’o 


The polytropie index y = 7/6 is determined by the stationary value of 6 =» 
7 . The entire range of observed 6 from Mariner 10 is 5.5 - 10.5 with 
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variability in the effective electron y of 1-10 - 1.22. This range of 
gammas is more nearly isothermal than adiabatic. This type of electron 
thermal behavior has been invoked previously by Burlaea et al. (1971)i 
Hundhausen and Montgomery (1971) and Hundhausen (1973). Recent theoretical 
modeling of two fluid solar wind dynamics by Goldstein and Jokipii (1977) 
and V. J. Pizzo (1978 private communication) using different polytrope 
indices for electrons and ions obtain: 1) better agreement with observa- 

tions than a one fluid approximation; and 2) tend to favor an electron y 
of approximately 1.15 - 1.2. This range of values contains the currently 
suggested value of y = 7/6 = 1.1 6 . Recently Feldman et al. (1978b) have 
suggested that solar wind electrons in stream interactions behave as if 
they obey a polytrope relation with y = 1.45, which is more nearly 
adiabatic than isothermal in character. We ‘will discuss these dynamical 
observations in detail in the last section of our paper, and suggest that 
their observations are consistent with nearly Isothermal electron behavior. 

The stationary value of 6 suggested by equation (4) implies a 
temperature variation when the density profile is specified. For large 
distances (n a r ), the condition of equation 4 (6 = 7) suggests that the 
thermal temperature varies as 

-1/3 

Tc(r)/Tc(ro) = (r/r^) , r > r^ (E* = e $) 


with a variability of exponent 0.-21 - 0.4 which results from the extremes 
of the observed 5. This range is similar to that suggested by Montgomery 
et al. (1968) based on Vela data taken at 1 AU and indirect arguments. 
Average radial variations of thermal electrons have recently been 
determined based on in situ observations at different radial positions. 
The reported power law fits in the approximate inverse square density 
regime are 

Mariner 10 

0.4 < r(AU) < 0.85 T -a r“°*^ +0.08 ogilvie and Scudder 

c 

(1978) 


Voyage r 2 

1.36 r(AU) < 2.25 T^ a Sittler ^ , 

(1978) 
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It is also interesting to note that this polytrope law implies an upper 
limit for the electron temperature of the corona as can be estimated along 
a given streamline as 


T„(1.03 R„) < T^(1 AU) 



(5a) 


The inequality arises due to the possible variation in y towards unity 

inside 10-20 R (cf. Parker (1963) for arguments). If we assume y = 7/6 to 
o 

R = 10 R and y = 1.1 between 1.03 - 10 R we can estimate 
o o 


n(1.03 Ro) 


T (1.03R ) = T^(1 AU) 
c o c 




n(10 R^) 


n(215 R^) 




(5b) 


”(10 Ro) 

Using composite empirical density profiles (Munro and Jackson (1977) and 
Sittler and Olbert ( 1 977 ) ) for coronal holes and in situ 1 AU electron 
temperature values given in Table 1 we obtain (via equation (5b)) 


Tg(1.03 Ro, Ug^ 740) = 8.9 x 10 ^°K (1.8 x 10^°K) 


T (1.03 R , U„„ ^ 400) = 1.7 X 10^°K (2.3 x 10^°K) 

0 0 oW 


where the value in parentheses is the upper limit via y = 7/6 independent 
of radius (equation 5a) . The intermediate computed values in Table 1 
appear in parentheses. These estimates are certainly consistent with the 
coronal hole electron temperatures inferred from soft X-ray data by Maxson 
and Vaiana (1977) of 9 x 10^ - 3 x 10^°K, but favor the lower portion of 
this range and possibly suggest that the temperature may be lower in the 
source region of the high speed wind than those for the slower wind. 

It is appropriate to discuss empirical estimates of e$. Feldman et al. 
(1975) have determined a number of parameters which characterize 1 AU 
electron distribution functions and have concluded (p. 4181) that "...halo 
electrons with energy greater than some breakpoint energy of >/' 60 eV" 

may be unbound. Comparable estimates of $ were obtained from the electron 
energy equation. This is to say that Z e$ (1 AU) is approximately -60 eV. 



20 


The quantity used to determine this number is highly quantized (of,., p. 
4183) with the effect that (a) Eg and its relative Eg^^ are systematically 
lower than a continuum definition would determine it to be and, (b) 
quantized in the proximity of 60 eV in increments of size 20 eV. The 
average thermal electron internal energy {per particle) is the average of a 
continuous parameter (Table 1 same paper) and is approximately 10,. 6 eV, 
which gives an empirical estimate of 


5 


Experimental 


+ 20 
60 ~ 0 


10.6 


eV 

= 6.6 ± 0.9 


which is clearly consistent with the theoretical value of 7. 

We thus conclude that relevant experimental estimates and theoretical 
arguments are consistent with the theoretical ansatz that ~ 7 kT^ 

derived from the properties of Coulomb collisions is synonymous with the 
polarization potential energy with respect to infinity provided 
equation (2b) is a valid approximation to the generalized Ohm's law along a 
tube of force from the observer to infinity. This is stated algebraically 
by 


®‘'Theory<''o> = 


( 6 ) 


Within this framework the polarization potential is the specific enthalpy: 


e$(r) = JL kT (r) = 7 kT^(r); y = 7/6. 

(Y-1) 

3. QRSERVATTOWAI. SUPPORT OF THEORETICALLY EXPECTED CORRELATIONS: STEADY 

STATE 


With the content of the transition energy in steady state clarified, 
the arguments for correlated behavior of different subpopulations of 
electrons with other solar wind variables are considerably simplified. We 
restrict our discussion in this section to steady state conditions. It 
should be clear from the contents of our discussion that the term 
"steady-state" should not be construed as meaning only the time-independent 
observational data gathered during some quiet time interval in the refer- 
ence frame of the spacecraft. Rather, the term implies the 
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time-independence as judged by an observer in a frame oo-rotating with the 
sun (or, more precisely, the frame in which the magnetic tubes of force may 
be considered as stationary if averaged over the fluctuations due to the 
possible presence of various wave modes) . Since the interplanetary medium 
is highly structured in heliocentric azimuth, all experimental data change 
in time, and only a careful, preparatory analysis can decide whether a 
given set of observations is steady in the sense defined above. If a data 
set can be declared as steady in the co-rotating frame and yet shows 
temporal variations in the spacecraft frame, a great variety of predictions 
can be obtained from a sequence of results of paper I, each referring to a 
different tube of magnetic field lines, characterized by a different state 
of plasma parameters within that tube. One can infer, for example, how 
observed changes in the magnitude of the wind velocity should be related 
(correlated or anti-correlated) to the theoretically expected changes in 
the various quantities characterizing the electron distribution. Such 
predicted correlations can be compared with the observed ones. Procedures 
of this type offer a strong test of the validity of our model. 

Correlations for Thermal Electrons 

As theoretically shown in the first section of this paper the sub- and 
lower energy transthermal electrons are in the best collisional contact 
with the ions and other electrons. The electrostatic quasi-neutrality 
condition also strongly couples the electrons and ions. It is therefore 
not surprising that the thermal electron density should emulate that of the 
ions and that the electrons should have a nearly common bulk speed with 
that of the ions. These phenomena have been reported previously starting 
with the work of Montgomery et al. (1968). Of course, the energy exchange 
rate between thermal electrons and protons is still very small so that 
their temperatures may be different. 

Whenever the local magnetic intensity is increased as a result of 
pressure imbalances across a flux tube, the transverse plasma pressure will 
be enhanced to the extent that energy losses along B allow. In the case of 
ions this loss is almost negligible and the heating proceeds essentially 
according to the double adiabatic hypothesis. By contrast, even subthermal 
electrons have sufficient speeds to remove energy from the compression 
regions . This mobility is counteracted by the changes in the local 
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polarization potential so as to maintain quasi-neutrality. These electrons 
nevertheless can conduct a large -amount of the compressive work away from 
the site of the compression. This feature of electrons was discussed from a 
local fluid standpoint by Hundhausen and Montgomery (1971) and Hundhausen 
(1973). The predicted polytrope behavior required in section 2, also 
allows a) a local description of thermal electrons and b) an equation of 
state that should imitate nearly isothermal behavior. 

Subthermal electrons do locally Coulomb scatter off of the embedded 
ions and other shielding electrons which on average are moving with the 
local fluid velocity. Therefore, if there are velocity gradient components 
(and therefore corresponding density and temperature variations) parallel 
to the local magnetic field with scale lengths comparable to, or longer 
than the electron mean free path, a subthermal electron can gain or lose 
energy depending on the sign of the gradient along the particle’s parallel 
velocity. In stream -regions this effect plus thermal conduction will tend 
to heat (if at all) the subthermal electrons in stream compressions more 
isotropically than one obtains by simple adiabatic compression. 

It is important to reemphasize that the work (+ or -) done on the 
parallel motions of the electrons can only be tapped by those electrons 
which recognize the existence of scattering sites moving with different 
bulk speeds. In this connection the subthermal electrons have a large 
probability of recognizing flow gradients along B with scale lengths L of 
the order of the thermal free path (1/4 - 1/2 AU at 1 AU) which is a 
typical scale length along the interaction pressure ridge at overtaking 
streams near 1 AU (V. J. Pizzo, private communication, 1978); by way of 
contrast the extrathermals have a 3o probability of Coulomb collisionally 
ignoring flow gradients on such a small scale compared to their free path. 

From this conclusion it appears necessary that the scale length, L, of 

/** 

bulk velocity gradients along the local tubes of force L = U*B/V(U’B) must 

be no shorter than x - in order that the electrons recognize the velocity 

mi p 

gradient and compress in phase with the ions; a magneto-fluid dynamic 
description with polytrope thermal electrons is only consistent in this 
circumstance. If this condition is violated (perhaps at shock boundaries 
where U and B are parallel throughout) then kinetic processes involving 
space charge, and time varying electric fields must be explicitly 
considered for such processes. 
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The radial variation of electron core temperature and electron bulk 

speed (among other parameters) as determined from the Mariner 10 PSE data 

are illustrated in Figure .6. These data span the radial range 0.45 - 0.85 

AU. The variation of the core temperature with radius has superposed 

structure caused by the dynamical interaction of corotating streams. This 

compressive effect expected of a classical gas can be clearly seen on the 

leading edges of the streams shown in Figure 6. If we concentrate on the 

trailing edges of these streams, where dynamical effects are minimal, there 

is a tendency for T^ to be lower locally when U is higher. This effect is 

consistent with the weak, inverse correlations between T and the bulk 

core 

speed reported by Serbu (1972) and Feldman ^ (1975) on the basis of 

1 AU observations. All of these correlations are consistent with expansive 
collisional cooling and compressive heating of a "classical gas". A 
complication for the former interpretation is that the electron tempera- 
tures in the coronal source regions of high speed flows .may be a priori 
lower than those for low speed flows, with the observations in the 
rarefactions being a mix of expansive cooling and initial (boundary) 
conditions. 

Correlations for Extrathermal Electrons 

The global, extrathermal component of solar wind electrons was the 
subject of several predictions in I. We suggested there that the local 
•bulk speed should be anti-correlated within a given corotating stream with 
fractional extrathermal phase density, extrathermal "differential 
temperature" and net heat flux. These correlations were not. local cause 
and effects, but were controlled by the electrical potential and the 
changes in collisional depth with energy of the observer relative to the 
distant collisional origin of the observed suprathermal phase density. The 
induced, local correlations given in the introduction were all consecuences 
of the global physics of the interplanetary polarization potential and. the 
integrated effects of the Coulomb cross section. We also suggested a 
theoretical reason that the halo temperature should be nearly independent 
of radius in the inner heliosphere In this section we organize the 
experimental support for these theoretical suggestions. 

In Figure 6 we also show the heliocentric variation of the "halo" , 
temperatures, the core Knudsen number, = VH. We call attention' to 
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the strong anti-correlation of the "halo'' differential temperature, 
with bulk speed, where 


, dln(f ) . 

= <-l [ Lr’>, 


E > E 


k - 


dE 


Data 


where the average < > is that implied by a linear least squares fit of Inf^ 
vs. E for E > E* in the proper frame. 

When the self-consistently determined bulk speed is locally high, is 
generally low and when the speed is lower, is generally larger. The 
principal isolated exceptions such as those at 0.634, and 0.485 AU are 
signals of stream interfaces discussed by Ogilvie and Scudder (1978) and 
can easily be screened out on the basis of and being locally 
anti-correlated and the large local core Knudsen number. Because of the 
radial variations superposed in Figure (6), it is best to rephrase our 
expectation that (within the stream) 


Sign (^Ih) = - Sign (^l^i_) 


dt 


dt 


in order that there be no confusion. There are at least three ways to 
understand Tj^ being smaller in high relative to low speed solar wind: 

( 1 ) perhaps there exists an omnipresent scattering agent , such as magnetic 
irregularities embedded in the accelerating flow; the suprathermals give up 
more energy (are "cooled") in scattering off these Irregularities in flows 
with larger spatial speed gradients (asymptotic bulk speeds higher, cf., 
e.g,, figure 6.1, Parker (1963)); or (2) extrathermal electrons at 1 AD in 
high speed flows have been most recently in collisional auasi-eauilibrium 
with cooler coronal plasma than is the case in low speed flows. (The site 
of collisional quasi-equilibrium of the 1 AU extrathermals near the corona 
may be (a) at different solar al'titudes with the underlying temperature and 
density profiles the same in different speed regimes as discussed in paper 
I, or (b) at the same coronal level with the electron temperature profile 
different (cooler) on high speed flux tubes than low ones.); or (3) supra- 
thermal electrons are expansively "cooled" by global Coulomb scattering at 
distances separated by sufficient radial distance that the local bulk 
speeds at the extremes of the Coulomb suprathermal mean free path are 
different; in this case the suprathermals give up energy to the flow via 
the embedded Coulomb scatterers . 
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From the measurements of whistler turbulence (Neubauer et al . (1977) 
and F. M. Neubauer, private communication (1978)) it is unlikely that there 
exists enough, appropriate turbulence to locally, “Fermi” cool the extra- 
thermal electrons in the case of alternative (1) above. The calculations 
of paper I, show that alternative (2a) is given quantitative support 
although (2b) cannot be completely ruled out. Alternative (3) is present 
in either (2a, b). Preliminary estimates suggest that (3) together with the 
hypothesis of coronal "source” regions of equal temperatures for high and 
low speed vfinds (2b) do not sufficiently discriminate between flow states 
to provide the observed cooling. We therefore conclude that (2a) and (3) 
above are consistent with theory and observations at the present time. 

In Figure 6 one can also see the theoretically suggested independence 
of (on average) with radius which has been .reported previously by 
Ogilvie and Scudder (1978) who reported a power law exponent consistent 
with zero. 

In Figure 7 we illustrate hourly averaged Mariner 10 PSE determinations 
of the halo/core density fraction and the bulk speed versus heliocentric 
radius between 0.45 and 0.9 AU. The theoretically expected 
anti-correlation of the extrathermal fractional density with local bulk 
speed is clearly shown. These data indicate that the halo energy regime is 
depopulated relative to the core in high speed flows when compared to low 
speed flows. We understand this as resulting from the higher potential 
energy barrier which the extrathermals must overcome in the high speed flow 
relative to the low speed flows. Since a strong , pervasive high frequency 
(ft) < whistler branch) turbulence field necessary to efficiently scatter 
halo electrons has not been directly observed in the solar wind (although 
surveys have been made, e.g. , Neubauer ^ al. , 1977), and as the effects 
discussed in paper I are unavoidable, it seems the explanation (2a, b above) 
of the extrathermal "cooling" and "depopulation" given in that paper is 
currently to be preferred. 

We consider in Figure 8 a superposed epoch analysis of electron 
properties through some "simple" very high speed streams observed at 1 AU 
as reported by Feldman ^ (1978a). These independent data also show 

the anti-correlation of Tj^ with the bulk speed. Of all the quantities 
plotted on this figure the halo fraction of the total density is the 

most volatile. We encourage the reader to inspect this record to deduce 



26 


the trend he perceives in this record for the stream behind the stream 
maximum. We have suggested by the solid line the trend we perceive in this 
data. The trend is certainly not positively correlated with the bulk 
speed. We suggest that this data also indicates an anti-correlation with 
speed with the halo fraction more populous in the lower speed region U 
400 than in the higher region. This is the same sense and magnitude of the 
(anti-) correlation inferred from the Mariner 10 data. 

In the. bottom panel of Figure 8 we call attention to the support these 
data provide for the expectation that the local heat flux should be in- 
versely related to the speed. The statistical study of Feldman ^ .al. 
(1976) at 1 AU support the generality of this situation. In this same 
connection Rosenbauer et al. (1977) also report that the heat flux is lower 
on average in high speed flows than in low speed flows at all radial 
distances sampled by the Helios spacecraft (0.3 < R ,(AU) < 1). 

We thus conclude that the observational data from three independent 
experimental groups on three different satellites which are relevant to our 
steady state calculations do provide strong support for the theoretical 
predictions made by us in paper I. 

4. OBSERVATIONAL SUPPORT OF THEORETICALLY EXPECTED CORRELATIONS: 

DYNAMICS 

Many of the changes seen in interplanetary electrons by an observer at 
fixed radial distance are associated with stream-stream interaction 
pressure ridges which are not necessarily corotating or, therefore, 

"steady” in the sense of the discussion of the previous section. We shall 
show in this section that the steady state theory developed in paper I, 
which has been extended and supported in the previous sections, should 
apply within dynamical stream interaction fronts. As an example of the 
theory's utility and self-consistency we shall illustrate, by direct 
calculation , the behavior expected of the extrathermal phase density at 
compressions associated with stream-stream interactions. The theory 
suggests that 1) the transition energy E*/kT is insensitive to dynamics 
and 2) that the extrathermal density variation at a proton compression 
depends sensitively on the polytrope index, y, appropriate for thermal 
electrons. The observations of the halo density variation at these ridges 
imply that y should have a value more nearly appropriate for isothermal 
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(jcl) as opposed to adiabatic- (5/3) behavior. On the basis of such 
successful application of the theory in dynamical situations, it is 
suggested that the theory for thermal solar wind electrons in the steady 
state contains the essential aspects of electron behavior even in the 
presence of relatively slow time variations represented by hydrodynamic 
disturbances . 

Superimposed upon the global steady state of the interplanetary medium 
(in the corotating frame) there exists some intrinsic evolution of the 
stream patterns , for example in response to boundary condition changes at 
the sun. The evolutionary re-adjustment of flow lines will result in 
temporal compressions and rarefactions of the steady state arrangement of 
material. Such “waves” will propagate with a speed less than or comparable 
to that of the speed of the wind, Uj the scale length of the stream-stream 
interaction fronts along B are typically L 1/4 - 1/2 AU at 1 AU (V. J. 
Pizzo, private communication, 1978). The intrinsic time scale t, of the 
variations experienced by the ambient plasma is of order L/U. This time 
scale is long when compared to the transit time of most electrons through 
the disturbance ; in other words such transient changes are slow as far as 
electrons are concerned, and the conclusions drawn in I from steady state 
considerations still apply in a semi-quantitative sense. Comparison of 
theory with observations of such transients may, therefore, be of further 
heuristic value. With the generalizations just made almost any stream 
interaction is suitable for comparison with the results to be calculated 
below — whether it is truly corotational or evolutionary. To illustrate 
this approach we address the problem of the response’ of the extrathermal 
density and transition energy E» at stream-stream interaction regions. 

Extrathermal “Compressions" 

Our approach is to note that within^compression (rarefaction) regions 
which have typical scale lengths along B larger than the local thermal 
electron mean free path that the subthermal (or core electrons) will behave 
as a classical gas — they will heat (cool) in phase with the increase 
(decrease) of ambient density to the extent consonant with energy 
transport. We have previously shown by steady state considerations that 
thermal electrons obey a polytropic law with a y = 7/6. We will specify 
the polytropic form but have y determined from experimental considerations. 
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where’ and are unperturbed reference values of P and n^, respectively, 
and the effective y is somewhere between the isothermal value of 1 and the 
adiabatic value of 5/3. 

At the same time, we note that from a collisional point of view that 
the halo electrons have a high probability of no.t collisionally recognizing 
the change in the density of thermal electron and ion scattering centers 
within a compression zone of scale length L ^ along the local tube of 
force. These particles will respond however to the perturbation in the 
polarization potential caused by the compression of the thermal electrons. 
The local change in the polarization potential energy of an electron 
relative to the ambient value at a nearby point R outside of the 
compression is (by generalized Ohm’s law in the spirit of approximations 
made in deriving equation 2b): 


eA$ = e($(r) -<J)(R)) 



n(r) 

In ( ), for (y = 1) 


k [T’^(r)-T^], for y 


( 6 ) 


where n(r) = n (R) + 6n(r) and T' (r) = T (R) + 6T(r); n , and T are 
assumed to be the spatially independent, nearby, unperturbed values outside 
of the compression along the given tube of force. Capital R denotes a 
place along the local tube of force removed from the compression but in the 
irnTncdiatp vicinity of it. 

We thus see that in a compression zone the polarization potential is 
larger than outside of the compression. This corresponds to a slight 
excess of positive charge in the zone. This additional attractive force 
accelerates the halo particles (which by the results of section 1 are to 
within 3 o) unimpeded by collisions to be observed at the kinetic energy 
E(r) = ^old^^^ + e[$ (r)-$ (R)] in the compressed region if they were at 
energy E^^^(R) outside of the compression region. 

We have schematically illustrated the electron distribution function 
prior to encountering a compression region in Figure 9(a) and within the 
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compression region in Figure 9(b). A Maxwellian in this format of Inf^ vs. 
E is a straight line. The differential temperature is proportional to the 
negative reciprocal of the slope of a straight line on this 
semi-logarithmic plot. The transition energy E* is operationally defined 
as the local proper frame kinetic energy where the core and halo Maxwellian 
parameterizations intersect. In general E*(r) will change in a systematic 
manner through the compression zone under the combined responses of the 
core's polytropic compression and the polarization potential changes that 
the halo experience. 

Since the halo "density" is equal to a partial moment of the 
distribution function: 


Hh = j di! I 
4ir 


f^ (w,ii)w^dw; 


(7) 


w* 


where w* is the speed determined by the transition energy, it is not clear 

how the halo "density" will vary with space Cor time for the observer r. - 

r + U (t-t )] through the compression. 

— o ~sw 0 

As shown in Figure (9a) the natural logarithm of the electron 
distribution function resembles 


fe>oor-e = '‘o ' ^ 

*'e>halc = \ - E/0=T^»)), (E > E«) 


(S) 


outside the compression, where A^,A^ are given by 

A„ = In (n„(R)/(TT^^^ wj)) 


A^ = In (nj^(R)/(ir^^^ 


where 1/2 m w ^ = kT , and 1/2 m w. ^ = kT. . From our physical arguments 
e c c 6 n 

regarding collisions and potential effects the modeled distribution 
function within the compression region at r is 
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^1" ^'e'W = “ E/{kT^), {E < E*’) 

(In f^*) = \ “ (E -eA$)/(kT^(R)), (E^E*') (9) 

halo 

where we have kept the halo temperature constant, since the dynamical 
compression will not , by our hypothesis , modify the "temperature” which the 
coronal conditions "impress" on the flux tube threading the compression 
region. For the purposes of these illustrative calculations we assume that 
the coronal conditions impress a constant in this region. [In principle 
the observed profile across adjacent sampled stream tubes could be sub- 
stituted to refine the calculation presented below. 3 The eA$(r) energy 
shift in the halo sub-population equation (9) models the local gain of 
kinetic energy which incident halo particles experience getting from R. to 
the observer at n. The new constant in equation (9) is given by 


= A,, -Mn [n^(r)/n^(Vr^’)3''^l 


where 


"^c/^ = (n^(r)/nQ)^“^ 

from the assumed polytropic law for the thermal electrons. 

To calculate the halo "density" variation through a model compression 

I 

(equation (7)) we must determine E* as a function of position along the 
tube of force. By solving (8) we obtain 

E«’(R) = BkT^, = T^(R) 

where T^ = T„(R) and 
o c 


3 = In [n /n • (T /T )^^^3 

T^-T 
h o 


using T^/T^ = 6, = 7 X 10^°K, and n^/n^^ = 19 at 1 AU we find 3 = 6.75. 

By solving (9) we obtain 
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E* (r) = 3kT» (r) [1 + 1/3 ln(n/n ], y = 1 

, 6T (5-3y) T. - 2yT 

E» (r) « 3kT* (r) [1 + (1 + )]. Y / 1. (11) 

(T^-Tq) 23(y-1)T^ 

It is important to note that the corrections terras to unity in the above 
two expressions are very small at 1 AU for typical excursions of the 
density and temperature. Thus, even in stream-stream interactions this 
model has the feature that the locally scaled break-point energy is only a 
weak function of slow dynamical effects. This is an important consequence 
of our model. It shows that dynamical effects do not impact the kinetic 
nature of the physics of E* which we introduced in the first section of 
this paper. It is therefore not surprising that a relatively narrow 
histogram, of E*/(kT^) such as obtained in Figure 4, could be obtained for 
a diverse collection of data in and out of evolving stream patterns. 

The variation of plasma quantities, scaled with distance i/.^ along the 
compressed flux tube, centered at maximum compression, according to this 
model are shown in Figure 10; [a is the arc length along the local tube of 
force and L the longitudinal scale length of the core compression]. Panel 
A shows the total density enhancement and the thermal electron temperature 
variation with arc length relative to their locally unperturbed values. 

The core density variation of panel D determines three variations of the 
thermal temperature in panel A referring respectively to the isothermal 
(y = 1), Y = 7/6 and adiabatic (y = 5/3) polytropic indices. The density 
and temperature profiles of panel (A) determine, via the generalized Ohm's 
law, the variation of an electron's potential energy along the tube of 
force under the three assumptions of the thermal polytropic indices. The 
extrathermals locally experience no resistance (collisions) and they will 
be attracted toward the compression and gain kinetic energy en route from 
the uncompressed local region of the flux tube. The core population does 
not have the mobility to overcome their shielding responsibilities of 
nearby ions. Only the extrathermals are clearly free to be attracted to 
the free charge within the compression. 

In panel C we show the computed variation of i^(i’) according to our 
model which entails evaluating (7) using equations (8-11). The extra- 
thermal response at maximum core compression of 200^ ranges from a 5%' 



32 


enhancement of the halo for near isothermal core behavior (a = 1, 7/6) to 
an l8/5 reduction for adiabatic core behavior (y = 1.45 yields a «/> 5^ 

of the halo) . In performing these calculations we have assumed 

5o 

that the halo had a constant temperature of 7 x 10 K and^that the 

5o 

unperturbed core temperature was 1.16 x 10 K. It is important to re- 
emphasize that the depicted variations are those an observer would report 
moving along the local tube of force through the compression. 

Before discussing measurements of the halo variations through stream- 
stream interaction regions , we should comment on the confusing and 
uncritical use in the literature of the word "compression", for any density 
enhancement which precedes the high speed flow immediately before the fast 
velocity increase observed in typical time records. (We do not mean the 
isolated non-compressive density enhancements structure in nearly uniform 
flow fields, discussed by Gosling ji,al.,, (1977).) A change in the number 
density, dn, may arise either from variations in size of a small volume 
element, d(AV), or changes in the number , (AN), of particles in that 
volume; 


AN d(AV) + d(AN) 
(AV)^ AV 


( 12 ) 


The variations with volume change are the compressive contributions to dn. 
The remaining contribution represent net change of number of particles in 
the volumes . As such this latter contribution to density variations 
represents the boundary conditions on the reference volume, (the halo 
density change computed above represents a combination of AV and AN). 

There is a growing body of evidence (experimental; Rosenbauer et al.. 1977 
and Gosling et al.. 1977, and theoretical; V. J. Pizzo, 1978, private 
communication) that a substantial portion of the density variation 
encountered prior to the steep stream velocity gradient is a remnant of the 
coronal boundary conditions in the initially denser, preceding slow stream 
which has been eroded by dynamical consequences of filling in the 
rarefaction of the preceding stream. 


The compressive zones are therefore not necessarily at the maximum 
density before the rise in the stream profile. Given that proton heat 



33 


conduction is almost negligible it seems warranted to locate the com- 
pression zone (where tubes of force are constricted) by the detection of 
simultaneously enhanced proton temperatures centered on the rapidly 
increasing slope of the bulk speed profile. 

Another intrinsic difficulty with examining temporal data records is 
that temporal variations do not translate into a unique picture of the 
spatial shapes of structures. This means, for example, that a density 
enhancement traversing the observer in one day at^a fixed solar wind 
evlocity H and imbedded magnetic field direction B, implies relatively 
little about the spatial extent of the structure along the local tubes of 
force . The density enhancement could, have almost any scale length along 6 
when compared with the apparent transverse scale length as determined from 
the width of structures convected over the observer. 

With these comments, we consider some published data from IMP-6 of the 
variations of the halo density through a stream- interaction region reported 
by Feldman (1975) (Figure 12) and reproduced here as Figure 11. The 
subpanels from top to bottom indicate the temporal records of the core and 
halo temperatures, the halo-thermal anisotropy, the core "drift” speed, the 
halo density and the break point energy Eg as defined in section 1 of this 
paper. The shaded zones denote from left to right, "...the region of large 
density compression where the added line indicates the position of the 
density maximum, and the high velocity region beginning with the approxi- 
mate position of the velocity maximum" These authors draw the following 
conclusions about the halo density data (p. 4191): "5. The halo density 
does not share in the core density compression which precedes the high 
speed region. Instead it is observed to be enhanced immediately after the 
main proton compression region. In this sense the halo behaves more like 
an independent, ubiquitous collisionless electron gas which responds to 
stream interaction structures induced in the solar wind plasma (similar to 
the response of low energy solar protons to interplanetary disturbances)." 
We have indicated the zone of maximum compression on the basis of proton 
temperature, Tp, peak with a solid dark box- above the halo density and core 
temperature traces. The halo density precipitously increases in the com- 
pression region we have determined using T . This zone is also in a region 
of enhanced but not maximal total density and near the maximum electron 
core temperature. Also note that Tj^ is nearly a constant across this 
enhanced halo density region. 
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From our model results the density of the halo can only be enhanced for 
nearly isothermal compressions y “ 1. Y - 7/6 also yields a halo enhance- 
ment; Y = 1.45 yields a reduction. This behavior appears to be at variance 
with the core temperature enhancement with enhanced density of the temporal 
records of Figure 1 1 which can be translated into a cross-sectional cut of 
the three dimensional variations at the stream interaction. However, as 
discussed previously, the cross sectional profile of a variable across 
different tubes of force does not determine the corresponding profile along 

/N 

B without additional assumptions. A nearly isothermal core electron 
profile along flow tubes in the crest of the compression (black zone) is 
consistent with the temporal record and our model. If we pursue -this line 
of reasoning, the heat transport along § from the compression away from the 
sun should be enhanced here in order to maintain nearly the same total 
electron temperature along the flux tube. This effect is not shown in this 
figure. However, this expected effect is clearly evident in Figure 8 
where, centered on a different halo density enhancement in the compression 
zone (defined by the protons black shading) , there is an unusually large 
heat flow. 

We thus conclude that the available data are consistent with the 
implications of our model of solar wind electrons even in the presence of 
slow variations . The steady state y which we used to interpret the 
polarization potential and suggest the distant T (r) profile also describes 
the slowly varying dynamic regime near compression regions. A portion of 
the T^ variation across adjacent flux tubes in the temporal records may 
therefore be interpreted as the remnants of coronal boundary conditions or 
the thermal history of the electrons en route to the point of observation 
as determined by the dynamics along those field lines which thread the 
compression pattern at some other radial distance. 

In view of our conclusion that thermal electrons obey a polytrope law 
with index 7/6, it is appropriate to test our results by examining a recent 
report by Feldman et al. (1978b) that concludes that Y £ 1.45 which is more 
nearly adiabatic (y = 5/3) than isothermal (y = 1). Their analysis was 
performed on the data partially displayed in our Figure 7 where these 
authors concentrated their attention on the fast compressed gas. They 
considered the observed temporal variation of the data which at 1 AU is 
principally across local tubes of forces. A crucial assumption of their 



35 


analysis is that the streamline constants at the footpoints of the sampled 
flux tubes have the same value on all sampled streamlines. The polytrope 
relation 


P(t,s)/PQ(0,s) =[p(t,s)/p(0,s)]^, (13a) 

where s indicates a given streamline, has been approximated as 

P(a,s) = 2^ Ep()1,s)]^; 2^ = constant (13b) 

where i.(s) denotes arc length along the streamline labeled s. The reported 
statistical analysis pertains to adjoining streamlines of compressed fast 
plasma, cf. figure 12(a) near the leading edge of the stream. Their 
analysis derives the effective least square which "characterizes” the 

data under the assumption that the variation of the streamline constant, 2, 
may be ignored on adjacent tubes of the force which thread the compression 
region and have been sampled. This assumption was made on the grounds that 
electron parameters are observed at 1 AU to be very homogeneous in the 
uncompressed high speed wind. 

As the thermal electrons dominate the partial pressure and density, the 

electron gas may behave in lowest order as a polytropic gas whenever the 

approximations leading to equation (2) are fulfilled (probably r > 20 R ). 

0 

If we take the empirical relation 

^eff 

Pg(t,s) a Pg(s) , = 1.45 ± 0.05 

we now inquire what variation of streamline constant S(s) is required to 
reconcile with y# = 7/6. If 


2(s) = P*/(p*)^,»p[j!,(s)] 


(Y*-Yeff^ 

f 


then these' author’s Ygff can be reconciled with the body of evidence 
compiled in this paper that y* “ 7/6 ±. .06. We thus determine that 
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Ya-Ygff = 7/6 - 1.45 = -0.28 + 0.08. 

Thus, for p[Jl(s)] varying by a factor of 5 as in these authors’ work, Z 
would have to have had a percentage variation of 


dz 2-s , 1 ■ » 

= </' (57 ± 8? (for z(s ) s > s ) 

« 

s z 

* 

where z is determined at the maximum compression within high speed flow 
» * 

(at s as indicated in Figure 12). Is it reasonable that Z be smallest on 

the most compressed flow tube? Since Z is a stream tube constant over a 

radial range, it seems appropriate to answer this question at the innermost 

boundary (R 20-30 R ) where we could defend by previous arguments that . 

w 

the same polytropic behavior obtains. The assumptions made by Feldman 
al . (1978b) is that the streamtube constants at the inner boundary are 
uniform and that z is a constant from heliocentric longitude i|)(s ) into the 
high speed wind (s > s*) even though these initial conditions have been 
impressed by the wind in which electrons are not behaving as a polytropic 
gas of the same y or may not be polytropic at all. There exists a body of 

g 

evidence which tends to suggest that pU = constant at 1 AU (5 ^ 1.5) 
(Burlaga and Ogilvie, 1970); mesa-like structures in flow velocity reported 
by (Rosenbauer ^ , 1977) at 0.3 AU show a similar anti-correlation of 

density, p, and solar wind speed, U. These mesa-like structures approxi- 
mate, but are not , square waves; the corners of the speed profiles are 
rounded as schematized in figure 12b. Assuming that = constant 
relation, it seems that the density profile will also have a complementary 
"shoe" to the descent from dense slow flow to sparse, fast flow. If we 
assume that the electron pressure is a constant across stream tubes, then 
we can relate precentage changes in Z to similar variations in speed as 

dZ ■ dU , 

3 > s P(s, 30 R ) = P (17) 

z U 

For the discrepancy between Tgff and y* to be attributed to flow tube 
constant variations , we require the density "shoe" to extend over such 
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at . 

longitudes (f)(s), s > s , from 63? of the mesa stream value until the speed 

profile reaches a flat maximum. This is a relatively thin and sharp speed 

stream. Under the alternate assumption that T^(s) is constant in longitude 

at 30 R , we require the density shoe to extend over the entire transition 
o 

from low to high speed flow; this is a very broad "stream-stream” boundary. 
The Helios data seem to indicate that high speed streams at 60 radii are 
narrow rather than as broad as the isothermal assumption requires. We have 
indicated in Figure 12b the percentage variations we require in order that 
values of y* = 7/6 and y^^^ may be reconciled with P (s, 30 R^) = P^. This 
is a plausible resolution of the interpretation of these data and our 
previous arguments in favor of y* = 7/6. 

6. DISCUSSION AND SUMMARY 

From a detailed consideration of Coulomb collision dynamics and the 
inhomogeneity of the solar wind, we have calculated the scaling of the 
transition kinetic energy, E*, of the local electron distribution for 
r > 0.45 AU to be 

E* - 7 kT^. 


We showed that this scaling is supported on average by all of the Mariner 
10 data (independent of radius) between 0.45 and 0.9 AU and also by a 
statistical analysis of published IMP (1 AU) data. Electrons with kinetic 
energies greater than E* are a special subpopulation of the suprathermal 
electrons that we have called the "extrathermals". The extrathermals have 
been shown to be the same group previously named the "halo" or "hot" 
component . 

By experimental and theoretical arguments we suggest that the 
transition energy, E*(r), and the polarization potential energy, e$(r) are 
on average synonymous in steady state for all radii greater than 0.45 AU 
where Hall effect and resistive contributions to $ are not important. A 
necessary condition that this result can be true is that the core electrons 
obey a polytrope condition 


P = 2(s)p^, y “ 7/6 + 0.06, 
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where S(s) are (perhaps different) streamline constants. An ancillary 
priviso for the use of the suggested polytrope behavior of thermal 
electrons as part of a locally neutral magneto-fluid-dynamic plasma 
description must also be understood? the scale length of the gradients of 
the component of the bulk velocity along the magnetic field should be long 
compared with the electron mean free path, 

When the density variation becomes inverse square with radius this 

polytrope law suggest that 

T =Ar"“, a = 1/3 
c 

which is consistent with recent in situ determinations on Mariner 10 (0.45 
- 0.9 AU) (Ogilvie and Scudder, 1978), Vela (Montgomerys^., 1968), and 
Voyager 2 (1.25 - 2.-3 AU) (Sittler S aj,. , 1978) of o = 0.3 + 0.08, 0.21 < 
a < 0.4, and a = 0.34 +. 0.l6, respectively. This asymptotic radial 
exponent of electron temperature variation is greater than, but very nearly 
2/7 as suggested by the fluid solution of Noble and Scarf (1963), Parker 
(1964) and less than the Whang and Chang (1965) solution which has an 
exponent of 2/5. This asymptotic variation also Implies that the mean free 
path for Coulomb scattering divided by the scale height along B, K| j , will 
become increasingly smaller and thus electrons in the ecliptic will become 
more fluid-like with increasing radial distance. If the same radial 
variations obtain for density and temperature over the solar poles, then 
K J I grows with radius as ^ indicating increasing departures from a 
Coulomb collisionally isotropized core and fluid behavior as r increases • 
over the pole. 

It is by no means clear, even now, whether the suggested equivalence of 

®e 

E»(r) with e$ (r) is a general result for a plasma with T^ » — T^ and an 

”p 

open magnetic topology when, resistive and Hall effects are negligible. We 
have shown that this condition implies a localized (polytrope) description 
for the thermal electrons which provide the majority of the shielding 
neutralizing charge density for the ions . In this sense this steady state 
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condition is nearly the requirement that a system of sufficient scale given 
sufficient time (steady state (t -»■ ®) ) will adjust itself so that the free 
charge is minimal. This association seems a posteriori reasonable in 
steady state since electrons with kinetic energy E > E are energetically 
unbound everywhere in the system and are correspondingly also allowed by 
Coulomb properties to respond to thermodynamic forces by being energeti- 
cally free with ballistic ranges against pitch angle changes long compared 
to the local mean free path. These particles, however, may be dynamically 
quasi-trapped via the cumulative backscattering they experience provided 
the system is large enough. The heliosphere appears to be large enough. 
This kinetic (local) and electrostatic potential (global) compatibility may 
be of more general importance than its application to the current problem. 

We have also shown theoretically that E*/kT^ is nearly independent of 
the slow dynamical evolution at stream boundaries. This feature of our 
model is an important consistency check of the fundamentally kinetic nature 
of the meaning of E* as derived solely from the properties of the Coulomb 
cross section. Had it been otherwise, we would have the unacceptable 
circumstance that the large scale dynamics of the solar wind affects local 
kinetic properties of the plasma, (Recall that E» was. defined in §1 as the 
energy below which escape from the local proximity was less than a 3a 
possibility.) 

A straightforward calculation with this model for conditions within a 
stream compression region near 1 AU shows that the extrathermal (halo) 
density can either increase or decrease within compression regions 
depending on the nature of the effective polytropic gamma of the local core 
compression. This sensitivity to the local nature of the core compression 
results from the interaction of the modified polarization electric field 
(determined primarily by the gradients of core T„, n ) with extrathermal 
electrons which are passing along the compressed tubes of force. In this 
sense the halo and core populations are locally coupled (even in slowly 
evolving stream interaction regions) through the strong force of the 
polarization potential . 

We concluded that, for two published examples of nj^^^o ^^^’i^tion at 
stream boundaries, the halo density probably does increase within stream 
interaction regions, contrary to the conclusions of Feldman et al. (1975). 
This conclusion is supported by our model calculations which predict halo 
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enhancements when the core compression is more nearly isothermal (e.g., 

7/6) than adiabatic {e.g. , 1 .45)~-together with the corollary that to 
maintain a nearly isothermal compression very large heat fluxes should be 
observed parallel to B within the compression region. This latter 
correlation is observed in the same paper. We have also suggested a way to 
reconcile their value with Y = 7/6 by modest streamline constant 
variations. There is thus no contradiction to the idea that the thermal 
electrons compress as if they were nearly isothermal (y = 7/6) at 
stream-stream interaction boundaries. 

In this paper we have presented strong observational support from data 
obtained on three different satellites reported by three independent 
experimental groups for all of the theoretically predicted correlations of 
paper I — namely : 

1) subthermal (E < kT^) electrons behave most nearly as a classical 
gas — such a gas compressively heats and expansively cools to the extent 
that heat transport permits and will participate in stream induced 
enhancements of the ions because of their strong collisional coupling to 
the ions and electrostatic plasma shielding requirements; 

2) the solar wind extrathermal fraction of the electron density is 
anti-correlated within steady state stream patterns with the local bulk 
speed (if.it is nearly the asymptotic value); 

3) the extrathermal electrons (with energy E > 7 kT^) do form a 
spectrally distinguishable subpopulation. The differential "temperature'* 
of these particles is observed to be anti-correlated with the local bulk 
speed. 

4) the heat flux carried by electrons' is anti-correlated with the bulk 
speed; and 

5) the extrathermal "temperature" is nearly independent of radius in 
the inner heliosphere. 

We thus conclude that the global and local Coulomb collisional effects 
discussed in paper I and elaborated here are essential aspects of the solar 
wind plasma as it is observed. The experimental support which the 
predictions of the theory in paper I enjoy, strongly suggests that the 
further development of that theory is warranted with the ultimate objective 
of understanding the impact that global heat transport has on the physics 
of coronal expansions. 
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FIGURE CAPTIONS 


Figure (1) Cuts of observed representative electron velocity 

distribution functions from Helios, fRosenbauer et al. . 
1977), ISEE-1 (Scudder and Olbert, 1978), and Voyager 
(Sittler ^ , 1978 ) showing the pervasive core-halo 

morphology of solar wind electrons. 


Figure (2) Thermal spread integral of "targets",^ (u) vs. scaled proper 
frame electron speed u = w/Wm lower horizontal axis, or 

1 j e 

u = ^/w_ upper horizontal axis. For the purpose of 
I- ,P 

illustration we have chosen w™ _/w™ =40 which corresponds 

i , e i , p 

to T ^ T . For all the considerations « (m„/m„) T„ 
e p p p e e 

1836 T , the structure of^ is as shown. 


Figure (3) Probability that a "test” electron will survive (right axis), 

[scatter (left axis)] while traversing the local root mean 

square target electron mean free path, as a function of 

scaled proper frame speed u = w/w.p^g. ^*xheory 

theoretical 3o definition of "survival". Other horizontal 

axes about this 3a level indicate 2a and 4a definitions of 

"survival". 2.19) and (= 2.21) are empirical 

determinations of the ‘ '.il- energy of f from IMP data 

and Mariner 10 (Mariner-Xe*ius-|ierGury, 1973) data — discussed 

below in text in detail. Stippled zone is the entire extent 

of the experimental deviations from the most probable values 

UviTTM* IMP data taken at 1 AU: MVM data between 0.45 and 0.9 
MVM 

AU. 


Figure (-4) Histogram of empirically determined breakpoint energy . 

E*„,„,(r) from Mariner 10 Plasma Science Experiment data 
MVM 

sealed by the local core temperature T (r); k is Boltzmann's 

c 


constant . 
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Figure (5) Correlation profile (data from Figure 18 of Feldman et al.. 

1975) of breakpoint energy Eg defined by Feldman et al. 

(1975) and core temperature. Optimal linear hypothesis of 
regression assuming errors in both variables superposed from 
the current analysis. 

Figure (6) ‘Radial variation of Mariner 10 electron data. Upper panel is 
electron ambient density estimates. Next panel contains core 
and halo temperature estimate. The bottom panel is a 
self-consistent estimate of the plasma bulk speed from 
electron data from the Plasma Science Experiment. 

Figure (7) Upper Panel: Variations of where asterisk denotes 

ambient estimates after spacecraft potential effects have 
been corrected. Lower Panel: Bulk speed variations from 

Mariner 10 Plasma Science Experiment electron observations. 

Figure (8) Composite synthesis from Figures 6,7> and 8 of Feldman et al. 

(1978b), showing superposed epoch variations of selected 
electron parameters near a corotating stream interaction 
region. Black bar in .speed panel denotes proton temperature 
enhancement . 

Figure (9) Schematic of electron distribution function response to a 
polytropic compression of the core electrons: (a) 

unperturbed; (b) within compression. 

Figure (10) Electron parameters through modeled compression region. 

Horizontal axis is arc length, z, along tube of force divided 
by the scale, L, of the density compression of the core. See 
text. 

Figure (11) Reproduced data figure from Feldman et al. , 1975, showing 

electron variations in compression regions. Black bar added 
to indicate the present authors’ opinion of where compression 
takes place, (note different location than density maximum, 
which Feldman et al. (1975) term as the compression). 
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Figure (12) (a) Schematic of density, P, stream tube constant, 2, and 

speedj U, across stream tubes observed at 1 AO. The 
horizontal axis is heliocentric longitude of the observer, 
(b) Hypothetical variations at R = 30 R^ which reconcile 
Yeff = 1.45 with an actual y« = 7/6. The "shoes" discussed 
in the text refer to the gradual transitions or ramps 
between flow states and are to be contrasted with the dotted 
square waves. 
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1.03 R 
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1.03 R_ 
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2 X 10 
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215 R„ 


4.8 
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ABSTRACT 


We have extended the theoretical considerations of Scudder and Olbert 
(1979) (hereafter paper I) to show from the microscopic characteristics of 
the Coulomb cross section that there are three natural subpopulations for 
plasma electrons: the sub thermals with local kinetic energy E < kT^; the 

transthermals with kT^- < E < 7 kT and the extrathermals E > 7 We 

present experimental support from^three experimental groups on three 
different spacecraft in the interplanetary medium over a radial range for 
the five interrelations projected in paper I between- solar wind electron 
properties and changes in the interplanetary medium: 1) subthermals 

respond primarily to local changes (compressions and rarefactions) in 
stream dynamics; 2) the extrathermal fraction of the ambient electron 
density should be anti-correlated with the asymptotic bulk speed; 3) the 
extrathermal "temperature" should be anti-correlated with the local wind 
speed at 1 AU; 4) the heat flux carried by electrons should be 
anti-correlated with the local bulk speed; and 5) the extrathermal differ- 
ential "temperature!' should be nearly independent of radius within 1 AU. 

From first principles and the spatial inhomogeneity of the plasma we 
show that the velocity dependence of Coulomb collisions in the solar wind 
plasma produces a bifurcation in the solar wind electron distribution 
function at a transition energy, E*. This energy is theoretically shown to 
scale with the local thermal (J' core) temperature as E*(r) _ 7 kT^(r). 

This scaling is observationally supported over the radial range from 0.45 
to 0.9 AU by Mariner 10 data and IMP data acquired at 1 AU. The extra- 
thermals, defined on the basis of Coulomb collisions, is synonymous with 
the subpopulation previously labeled in the literature as the "halo" or 
"hot" component. If the transition energy should be required to equal the 
polarization potential energy, e$(r) , for all radii beyond the observer, 

(as motivated by Mariner 10 data) the thermal electrons should obey a 
polytrope law with index y = 7/6. In this circumstance the polarization 
potential is equal to the specific enthalpy: e$(r) = [y/(Y-l)] kT (r) . 

This relation probably does not obtain in the corona, inside r < 10-20 
(cf. Paper I). In the asymptotic spherically symmetric solar wind the 
inverse power law index for radial -variation required for the thermal 
electrons should be a = 1/3 which is consistent with the recent in situ 
determinations between 0.45 and ^f>3 AU. We thus provide the first 
self-consistent argument for associating E*(r) with e$(r) . This 
theoretical asymptotic (r -► ») thermal electron temperature variation is 
between the conduction dominated (2/7) and the inviscid solution (2/5). In 
the proximity of the ecliptic the ratio of the thermal electron Coulomb 
mean free path to scale length is shown to decrease with increasing radial 
distance. By contrast, over the solar poles the same asymptotic tempera- 
ture and density variations imply that the Coulomb thermal mean free path 
will increase with increasing radial distance. 

A straightforward extension of our model to slow time dependent 
situations shows that the fundamental scaling of E* with temperature and, 
by inference, the polytrope law for thermal electrons is theoretically 
maintained within compression and rarefaction regions. An illustrative 
calculation is given of the extrathermal density's response in stream 
compressions near 1 AU. These calculations show that the extrathermal 
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electrons density can either increase or decrease depending on whether the 
core electrons compress more nearly as an isothermal gas or an adiabatic 
one. On the basis of these calculations and a reexamination of published 
data we suggest that the extratherraal density is probably enhanced within 
stream compression regions. This implies (as observed) that (a) the core 
electrons should compress more nearly as an isothermal (y = 1) rather than 
as an adiabatic (y = 5/3) gas> and (b) a very large electron heat flux 
should be observed leaving the compression region along the magnetic field. 

We thus conclude that the global and local Coulomb effects discussed in 
paper I and elaborated here are essential aspects of the solar wind plasma 
as it is observed. The experimental support which the predictions of the 
theory in paper I enjoy, strongly suggests that the further development of 
that theory is warranted with the ultimate objective of understanding the 
impact that global heat transport has on the physics of coronal expansions. 



